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ABSTRACT 

We present a cross-correlation analysis of the WMAP cosmic microwave background (CMB) temper- 
ature anisotropies and the SDSS galaxy density fluctuations. We find significant detections of the 
angular CMB-galaxy correlation for both a flux limited galaxy sample (z ~ 0.3) and a high redshift 
(z ~ 0.5) color selected sample. The signal is compatible with that expected from the integrated 
Sachs- Wolfe (ISW) effect at large angles (8 > 4°) and the Sunyaev-Zeldovich (SZ) effect at small 
scales (9 < 1°). The detected correlation at low-z is in good agreement with a previous analysis using 
the APM survey [z ~ 0.15). The combined analysis of all 3 samples yields a total significance better 
than 3<7 for the ISW and about 2.7a for the SZ, with a Compton parameter y ~ 10 -6 . For a given flat 
ACDM model, the ISW effect depends both on the value of f^A and the galaxy bias b. To break this 
degeneracy, we estimate the bias using the ratio between the galaxy and mass auto-correlation func- 
tions in each sample. With our bias estimation, all samples consistently favor a best fit dark-energy 
dominated model: J7a ~ 0.8, with a 2a error Qa = 0.69 — 0.86. 
Subject headings: cosmic microwave background, cosmology: observations 
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1. INTRODUCTION 

A recent study ijFosalba fc Gaztanaga 1200*31) has cross- 
correlated the cosmic microwave background ( CMB) 
anisotropies measured by WMAP ijBennett et alJl2003|) 
with galaxy fluctuations in the APM Galaxy Sur- 
vey l)Maddox. Efstathiou. Sutherland, fc Lovedavlll990D 
to find significant detections for both the integrated 
Sachs- Wolfe (ISW) and thermal Sunyaev-Zeldovich (SZ) 
effects. The ISW detection is in agreement with 
other analyses based on X-ray and r adio sources 
(iBoughn fe Crittenden! I2003t iNolta et all 120031) w hile 
Hcrnandez-Monteagudo & Rubino-Martinez (2003) fail 
to detect the SZ effect when com paring WMAP to d if- 
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ferent cluster templates (see also 
It should be stressed nevertheless that cluster or galaxy 
group catalogues are too sparse and typically produce 
worse signal-to-noise ratios than galaxy surveys. More- 
over, depending on the sample there could be a sig- 
nificant cancellation of the ISW and SZ effects on 
scales smaller than a few degrees (see 2J. In this let- 
ter we cross-correlate the WMAP CMB temperature 
anisotropies wi th galaxies from the Sloan Digital Sky 
Survey (SDSS; I York et al .1 120001) . When we were fin- 
ishing this work we b ecame aware of a similar analysis 
(jScranton et alJl2003() that uses different color and pho- 
tometric redshift selected samples from the SDSS. 

2. DATA 

We make use of the largest datasets currently avail- 
able to study the CMB-galaxy cross-correlation. In or- 
der to probe the galaxy distribution, we have selected 
subsamples from the S PSS Data Release 1 (SDSS DR1; 
lAbazaiian et alJl2003|) which covers ~ 2000 deg 2 (i.e, 5 
% of the sky). The samples analyzed here have differ- 
ent redshift distributions and a large number of galaxies 
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(10 5 -10 6 , depending on the sample). We concentrate our 
analysis on the North sky (~ 1500 deg 2 , ie, 3.6 % of the 
sky), because it contains the largest and wider strips. 
The South SDSS DR1 (~ 500 deg 2 ) consists of 3 narrow 
and disjoint 2.5° strips, which are less adequate for our 
analysis. 

Our main sample, hereafter SDSS all, includes all ob- 
jects classified as galaxies with extinction corrected mag- 
nitude r < 21, and a low associated error (< 20%). This 
sample contains ~ 5 million galaxies distributed over the 
North sky. Its predicted redshift distribution is broad 
and has a median redshift z ~ 0.3. Our high-redshift 
sample (SDSS high-z thereafter) comprises ~ 3 x 10 5 
galaxies, with z ~ 0.5. It was selected by imposing mag- 
nitude cuts and color cuts perpendicular to the redshift 
evolution and the spectral type variations based on theo- 
retical spectral synthesis models. We shall also compare 
our re sults to the APM analysis in lFosalba fc Gazta naga 
(2003), who used a bj — 17 — 20 sample, z ~ 0.15, area 
~ 4300 deg 2 and 1.2 million galaxies. 

For the C MB data, we use th e first-year full-sky 
WMAP maps l)Bennett et al.ll2003D . Since the observed 
CMB-galaxy correlation is prac tically independent of 
the W MAP frequency band used l)Fosalba fc Gazta naga 
2003), we shall focus on the V-band (~ 61 GHz) as it con- 
veniently combines low pixel noise and high spatial res- 
olution, 21'. In addition, we have also us ed the W-band 
and a foreground "cleaned" WMAP map l|Tegmark et alJ 
2003) to check that our results arc free from galactic con- 
tamination. We mask out pixels using t he conservative 
KpO mask, that cuts out 21.4% of the sky l)Bennett et alJ 
2003). All the maps used have been digiti zed into 7' pix- 
els using HEALPix 4 , (jGorski et alJl!999|) . 

3. CROSS-CORRELATION AND STATISTICAL TESTS 

We follow the not ation introduced in 
iFosalba fc Gaztanaga I l)2003[) . We define the cross- 
correlation function as the expectation value of density 

4 http:/ /www. eso.org/science/healpix 
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Fig. 1. — Errors in the cross-correlation wtgW 1™ the dis- 
persion in 200 Monte-Carlo simulations (solid line) as compared 
with the mean and dispersion (squares with errobars) in the Jack- 
knife (JK) error estimation over the same simulations. Dashed line 
correspond to the JK error in the real WWLAP-SDSS all sample. 



fluctuations Sq — Nq/ < Nq > — 1 and temperature 
anisotropies At = T — Tq (in (jK) at two positions h\ 
and fi2 in the sky: wtg{9) = (Ar^ijfc^)), where 
6 = |n a - fii|. 

We compute the CMB-galaxy correlation and the as- 
sociated statistical er ror-bars using the jack-knife (JK) 
method described in ijFosalba k, Gaztanaea 1 12003^) and 
references therein. The survey is divided into M = 16 
(we find similar results for M = 8) separate regions on 
the sky, each of equal area. The wtg analysis is then 
performed M times, each time removing a different re- 
gion, the so-called JK subsamples. The covariance CV,- 
for wtg between scales 9i and 9j is obtained by re-scaling 
the covarian ce of the JK subsample s by a factor M — 1 
(see Eq.[3] in lFosalba fc Gaztanagal 1)20031) ). To test the 
JK errors and covariance we have also run 200 WMAP 
V-band Monte-Carlo (MC) realizations. We add random 
realizations of the measured WMAP tem perature angu- 
lar power-spectrum l)Bennett et alJl20 03) to those of the 
white noise estimate d for the relevant frequency band 
(Hins haw et al.ll2003f) . For each MC simulation we esti- 
mate the mean "accidental" correlation wtg of simulated 
CMB maps to the SDSS galaxy density fluctuation map. 
We also estimate the associated JK error in each MC sim- 
ulation. Fig^ compares the 'true' sampling error from 
the dispersion of wtg (9) in 200 MC simulations with the 
mean and dispersion of the JK errors over the same sim- 
ulations. The JK error gives an excellent estimate of the 
'true' error up 9 ~ 5 degrees. On larger scales it only 
underestimates the error by 10 — 20%, which is hardly 
significant given the uncertainties. 

Fig0 shows wtg for the different samples together 
with the corresponding JK error. It turns out that the 
the JK errors from the real WMAP sample are in some 
cases smaller (up to a factor of two) than the JK er- 
rors (or sample to sample dispersion) from the MC sim- 
ulations. Fig|2 shows, as a dashed-line, the comparison 
for the SDSS all sample, which exhibits the largest dis- 
crepancy. This difference in error estimation is not to- 
tally surprising as the MC simulations do not include 
any physical correlations and use a CMB power spec- 
trum that is valid for the whole sky, and not constraint 
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Fig. 2. — WMAP-SDSS correlation: long dashed-line shows the 
measurement for the SDSS HIGH-Z sample, while the solid line 
displays the correlation for the SDSS ALL sample. For reference, 
short-dashed line displays the same measurement using the APM 
galaxy survey instead of SDSS. Boxes show 1-cr error-bars. 



as to match the CMB power over the SDSS region. The 
JK errors provide a model free estimation that is only 
subject to moderate (20%) uncertainty, while MC errors 
depend crucially on the model assumptions that go into 
the simulations. Despite these differences in the MC er- 
ror estimation the overall significance for the detection 
turns out to be smiliar, as explained in §4.1. 

We derive the significance of the detected correlation 
taking into account the large (JK) covariance between 
neighboring (logarithmic) angular bins in survey sub- 
samples (but see also §4.1). Adjacent bins at large scales 
(9 > 4°) are correlated at the ~ 80% level, dropping to 
~ 40% for alternative bins. Bins at smaller scales are 
progressively more correlated. To assign a conservative 
significance for the detection (ie against wtg = 0) we es- 
timate the minimum x 2 At for a constant wtg and give 
the difference A% 2 to the wtg = null detection. For ex- 
ample, at scales 9 = 4—10° we find: wtg — 0.53±0.21^K 
for the SDSS high-z sample, wtg = 0.26±0.13/iK for the 
SDSS all sample and w T g = 0.35±0.13AiK for the APM, 
in all cases we give 1-tr errorbars. 

We find the largest significance in the CMB-galaxy cor- 
relation for the SDSS high-z sample: A% 2 =9.1 (ie prob- 
ability, P = 0.3% of no detection) for 9 < 10° (being 
Xmin = 14.6 for wtg = 0.55/iK with 11 d.o.f., although 
the fit is only approximate as the signal drops with scale) . 
In order to assess the significance levels for the ISW and 
SZ effects from the observed CMB-Galaxy correlations, 
we shall first introduce model predictions. 



4. COMPARISON TO PREDICTIONS 

The temperature of CMB photons is gravitationally 
rcdshifted as they travel through the time-evolving dark- 
matter gravitational potential wells along the line-of- 
sight, from t he last scattering su rface z s = 1089 to 
us, z = l)Sachs & Wolfd 11967ft . At a given sky 
position n: AT ISW {n) = -2 J dz $( n, z), and for 
a flat univers e V 2 $ = -AnGa 2 p m S (see Eq.[7.14] in 
Peebles (1980)). In Fourier space it reads, <&(k,z) — 
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Fig. 3. — Theoretical predictions: (Bottom panel) Continuous, 
long and short dashed lines show the ISW, SZ and lensing pre- 
dictions. Different sets of lines correspond to the APM (black), 
SDSS all (red) and SDSS high-z (blue) samples. (Top panel) The 
total prediction (ISW+SZ+Lensing) for the 3 samples. We have 
assumed a ACDM model with a fixed b ga s = 2 in all cases, 6 = 3 
for SDSS high-z and b = 1 for APM and SDSS all. 



-3/2n m (H Q /k) 2 S(k,z)/a, and thus: 

w T s G w (e) =< A^Sa >=J§ P(k) 9 {k8) (1) 

being, g{k6) = l/2n J dz W ISW {z) W G (z) j (k8r), 
where the ISW window is Wisw — —3£l m (H /c) 2 F(z), 
with c/Hq ~ 3000/i Mpc" 1 , F = d(D/a)/dr = 

(H/c)D(f — 1), and / ~ flm^iz) quantifies the time 
evolution of the gravitational potential. The galaxy win- 
dow function is Wq — b(z) D(z) (/>g(z), which depends 
on the galaxy bias, linear dark-matter growth and the 
galaxy selection function. The ISW predictions for the 3 
samples are shown in in bottom panel of FigEI Unless 
stated otherwise, we use the concordance ACDM model 
with n m = 0.3, n A = 0.7, T ~ M2 m = 0.2 and cr 8 = 1. 

The weak lensing effect prediction is quite simi- 
lar to the ISW, we just need to replace the time 
derivative of the Newtonian potenti al by its 2D 
Laplacian (|Seliak fc Zaldarriagal 11996ft: W Lens = 
3k 2 rt m (Ho/c) 2 (D/a)/d(r) d{r) being the angular dis- 
tance to the lensing sources (with: d(r s — r)/d{r s ) ~ 1). 

For the thermal Sunyaev-Zeldovich (SZ) effect, we as- 
sume that the gas pressure 8 gas fluctuations are traced 
by the galaxy fluctuations 5 gas ~ b gas Sg with a rela- 
tive amplitude given by the gas bias, b gas ~ 2, repre- 
sentative of galaxy clusters, although b gas is uncertain 
to within 50 % on linear s cales and for low-z sources 
(iRefr eeier fc Tevssierl20 021 . A si mple conservative est i- 
mate of the SZ effect is given by ijRefregier et alJ l2000): 

w tg(@) = ~b gas AT wgg{9) (2) 

where AT is the mean temperature change in CMB pho- 
tons Compton sca ttered by electr o ns in hot intraclus- 
ter gas. Following Re fregier et alJ 1)200(1) . we calculate 
AT = j{x)y T , where T ~ 2.73K is the mean CMB 
temperature, y is the mean Compton parameter induced 
by galaxy clusters, and j(x) = —4.94 is the negative SZ 



spectral factor for the V-band. The Compton parameter 
can be calculated integrating along the line of sight the 
normalized galaxy redshift distribution convolved with 
the volume-averaged density-weighted temperature. The 
latter is obtained from the mass function and the M-T 
relation. We assume the Seth &: Tormen m ass function 
(ISheth fc Tormen Ill999l ISheth et al ll20f)l and the M- 
T relation given bv lBorga^rTeT^L 1 l)1999|) In summary, 
we obtain for the WMAP V-band, AT = 6.65/xK for the 
SDSS all sample, and AT = 6.71^K for the SDSS high-z 
sample which corresponds to y ~ 1.35 x 10~ 6 for both 
samples. The SZ predictions for the 3 samples are shown 
in the bottom panel of Fig|3| Note the galaxy auto- 
correlation explains most of the differences observed. 

The total predicted correlation is thus the sum of three 
terms: the ISW, thermal SZ and Lensing contributions, 
wtg = w t S g W +WTG JrW TG aS ■ FigOJshows individual con- 
tributions of these effects (bottom panel) and the total 
(top) for the 3 samples analyzed. The ISW effect typi- 
cally dominates for angles 8 > 4°, while the SZ effect is 
expected to be significant on small scales 8 < 1°. Lensing 
is found to be negligible at all scales for our samples. 

Before we can make a direct comparison between the- 
ory and observations, we shall address the issue of galaxy 
bias. The higher redshift sample requires a high bias 
b > 1 to explain the large cross-correlation seen at all 
scales (the SZ effect being smaller at high redshift). At 
low redshifts the measured correlation is dominated by 
the thermal SZ on small scales (6 < 1°) and by ISW on 
large-scales (8 > 4°). Here no bias is required to match 
the observations. This agrees quite well with our self- 
consistent bias estimation: for each sample we can esti- 
mate the ratio b 2 ~ u>gg/ w mm , where wmm and wgg 
are the (theoretically predicted) matter and (measured) 
galaxy auto-correlation functions. For APM and SDSS 
all samples we find b 2 ~ 1, while for the SDSS high-z 
sample we get b 2 ~ 6. 

4.1. Significance tests 

ISW effect: On large scales 8 > 4°, the ISW effect is 
expected to dominate for all survey depths (see Fig[3J|. 
Therefore, from the large-angle CMB-Galaxy correlation, 
we can directly infer the ISW effect (ie, wtg — w t S g W > 
see end of <J3). In particular, for the SDSS high-z sample, 
a constant correlation fit rejects the null detection with 
high significance A^ 2 = 6.0 (P = 1.4%), comparable to 
the level found for the APM survey, A\ 2 =6.1 (P = 
1.3%). A smaller significance is obtained for the SDSS 
all sample: A\ 2 = 3.9 (P= 4.8%). Alternatively, we can 
use the uncorrelated MC simulations (see §3) to have an 
independent estimate of the significance. When a par- 
ticular MC simulation has an accidentally large value of 
wtg it also has a large JK error associated. We can thus 
assign a significance to our measurement by asking how 
many of the 200 MC simulations have a value of wtg 
equal or larger than the observations with an assosiated 
JK error equal or smaller than found for the observa- 
tions. We find that only two of the MC simulations fulfil 
this condition in any of the samples, meaning that the 
significance of each detection is better than 1% (for each 
of the 3 different data samples). 

Since these samples are basically independent, we can 
combine them to infer a total significance for the ISW 
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Fig. 4. — Estimating dark-energy: Long-dashed, short-dashed 
and dot-dashed lines show the probability distribution for Qa in 
the SDSS all, APM and SDSS high-z samples. The combined dis- 
tribution (for 3 d.o.f.) is shown by the solid line. 



detection: we find a total A X 2 = 16 (P = 0.1% for 3 
d.o.f) corresponding to a 3.3er. Note we could do bet- 
ter using a (scale dependent) ACDM model theory pre- 
diction, but at the cost of introducing model dependent 
detection levels. Moreover, we can further include the 
ISW-dominated small angle bins in our deepest sam- 
ple, where SZ is negligible, increasing the significance 
to Ax 2 = 18.8, (P = 0.03% for 3 d.of.), ie we detect the 
ISW effect at the a 3.6a level. 

SZ effect: We can estimate the significance of the drop 
in the signal at small angles in the SDSS all and APM 
samples due to the SZ effect (see Fi^3Jl using the best-fit 
constant at large angles (ie, the ISW signal) and ask for 
the observed deviation from such value at smaller scales. 

0.27 ± 0.11 for SDSS all, 
0.41 ± 0.16 for APM (1-cr errorbars). Note 
this is conservative because the ISW increases slightly 
as we approach smaller scales (see §IV). This test gives 
A X 2 = 5.5 (P = 2%) for the SDSS all sample and A X 2 = 
8.5 (P = 0.3%) for the APM. 



For 9 < 1°, we find: w§,q = 
and Wj>g 



5. DISCUSSION 

We have measured the CMB-galaxy correlation using 
WMAP and the SDSS DR1 galaxy survey. We measure 



a significant cross-correlation at low (z ~ 0.3) and high 
(z ~ 0.5) redshifts. We detect a positive correlation on 
large-scales induced by the ISW effect at the 2a level 
for the (broadly distributed) low-z sample. This corre- 
lation is similar to that m easured for the lower r edshift 
[z ~ 0.15) APM galaxies l|Fosalba fc Gaztanaga 1 82003) . 
although the latter has a larger significance, 2.5er. More- 
over, the significance of the detection raises to 3cr for the 
SDSS high-z sample. The combined analysis for the 3 
samples gives a 3.6<r significance (see H4.1|l . 

Our measurements at large scales are in good agree- 
ment with ISW predictions for a dark-energy dominated 
universe. Fig0]shows the probability distribution for fl\ 
in a flat ACDM model. We have fixed a s = 1, h = 0.7 
and Qjvf + = I- As we vary Q\ the shape parameter 
for the line ar power spectrum P(k) c onsistently changes 
T = hCt M l)Bond k. Efstathioulll98l . We only use the 
data for 9 > 4°, where the ISW is the dominant contri- 
bution. We fix the bias b by comparing the matter angu- 
lar auto-correlation function in each model to the galaxy 
auto-correlation in each sample. We find b ~ 1 for the 
APM and SDSS all samples and b ~ for the SDSS 
high-z sample. The A% 2 value in each model refers to the 
minimum \ 2 fit to a constant in the range 4° < 9 < 10°. 
As can be seen in the figure, all samples prefer large val- 
ues of ft a, with the best fit J7a ~ 0.8 with a 2a range 
A = 0.69- 0.87. 

We also see evidence (2.7 a level) for the thermal SZ 
effect from the drop of the CMB-galaxy correlation on 
small-scales in the low-z samples of SDSS and APM 
galaxies. These new measurements can be used to con- 
strain the redshift evolution of the physical properties of 
gas inside galaxy clusters. 
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